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Abstract 

Spatial qudit states can be realized by using multi-slits to discretize the transverse momentum 
of a photon. The merit of this kind of spatial qudit states is that the implementation of higher 
dimensional qudits is relatively easy. As we have recently shown, the quantum states of these 
spatial qudits can be analyzed by scanning a single interference pattern. This method of single 
scan tomography can also be applied at higher dimensions, but the reconstruction becomes more 
sensitive to smaller details of the scanned patterns as the dimensions increase. In this paper, we 
investigate the effect of finite measurement resolution on the single scan tomography of spatial 
qutrits. Realistic measurement operators describing the spatial resolution of the measurement are 
introduced and the corresponding pattern functions for quantum state reconstruction are derived. 
We use the pattern functions to analyze experimental results for entangled pairs of spatial qutrits 
generated by spontaneous parametric down-conversion (SPDC). It is shown that a reliable recon- 
struction of the quantum state can be achieved with finite measurement resolution if this limitation 
of the measurement is included in the pattern functions of single scan tomography. 
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I. INTRODUCTION 



In quantum information science, information processing tasks such as quantum cryptog- 
raphy, dense coding, and teleportation can not only be implemented using two-dimensional 
qubits, but are sometimes more efficiently performed using ci-dimensional qudits as carriers 



of information 



are sometimes more emcientiy periorn 

ion Q, y, y, Q, y, y, y, y, y, y 0, q. 



In optical implementations, we can make 



use of the degrees of freedom of photons to define the qubits or qudits. While photon polar- 
ization is a two-dimensional degree of freedom, the spatial (transverse momentum or angular 
momentum) and spectral (time-frequency) degrees of freedom are intrinsically continuous 
and can be used to define arbitrarily dimensional qudits by appropriate discretization. For 
example, qudits can be defined by selecting a set of d angular momentum eigenstates rep- 
resented by Gauss-Laguerre modes [a, l9||. Alternatively, multi-slits can be used to define d 
transverse modes of the light field 10| . In this case, all modes are essentially equivalent and 
an increase of d by adding more slits is relatively easy. Therefore, the method of using multi- 
slits to discretize the photons' transverse momentum is very promising for the realization of 
qudits with higher dimensions. 

One of the challenges in realizing higher dimensional qudit systems is that the higher 
dimensional states require more precise and detailed measurements for analysis and control. 
For this purpose, it is useful to adapt the method of quantum state tomography to the 
specific qudit system 
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Recently, we have presented a method of 



quantum state tomography for rf-slit spatial qudits based on the spatial patterns obtained 
by scanning the position of photon detection in a plane between the focal and image planes 



of a lens 



171 ]. In this initial work, we demonstrated the method for double-slit qubits and 



pointed out that the extension to higher dimensions is straightforward. However, as the 
following results for triple-slit qutrits show, an increase of the slit number d requires an 
increased spatial resolution to observe the theoretically predicted patterns. Since such high 
spatial resolution is difficult to obtain, it is better to include the finite resolution of photon 
position detection in the method of quantum state tomography. It is then possible to obtain 
reliable results even though the measurement setup cannot fully resolve the details of the 
interference patterns. 

In this paper, we present the formalism for single scan tomography with finite spatial 
resolution and demonstrate it experimentally for the case of spatial qutrits. For this pur- 
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pose, we describe the measurements in terms of realistic measurement operators including 
the resolution of the detection setup. The elements of these measurement operators then 
provide the actual pattern functions observed in the experiments. The quantum state is 
reconstructed by fitting the set of pattern functions to the detection distribution. The effect 
of finite resolution is clearly seen in the improvement of the fit achieved by using the realistic 
measurement operators instead of the theoretical infinite resolution interference patterns. 

The rest of the paper is organized as follows. In Sec. [Tll we review the principle of 
single scan tomography and show how finite resolution can be represented by the measure- 
ment operator. In Sec. IIIIl we apply the method to two triple-slit qutrits measured by 
photon detection in a plane between the focal and image planes of a lens. We derive the 
pattern functions with and without finite resolution and illustrate the effects of finite res- 
olution. In Sec. IIVI we describe our experimental setup and present measurement results 
that demonstrate the entanglement of the qutrit pairs. In Sec. |Vl we present and discuss 
the experimental results for single scan tomographies of spatial qutrit states prepared by 
conditional measurements on the other qutrit. The complete density matrix of the two 
qutrits is then reconstructed from a sufficiently large set of conditional measurements. Sec. 
IVll concludes the paper. 

II. SINGLE SCAN TOMOGRAPHY WITH FINITE SPATIAL RESOLUTION 

Multi-dimensional spatial qudits can be defined by discretizing photon transverse mo- 
mentum using multi-slits. If the quantum state of the photon passing through the z-th slit 
is given by \si), a (i-slit generates a qudit described by the orthogonal basis After 
the state of the photon is discretized by the slits, the photon propagates in continuous free 
space. The state of the spatial qudit can then be analyzed by measurements of photon 
position in this continuous space. Specifically, the probability of detection at a point x in 
any plane behind the d-slit is given by a measurement operator M{x) that acts on the d- 
dimensional Hilbert space of the spatial qudit. For a qudit state given by a density matrix 
p, the probability of detection is 




(1) 
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where Mij{x) and pji are the matrix elements of the operators in the basis of the spatial 
qudit Thus the detection probability P{x) is a linear combination of the measurement 

operator elements with the coefficients given by the density matrix elements of the qudit 
state. If all of the pattern functions Mjj(x) are linearly independent, it is possible to invert 
and solve the equation for the density matrix elements, so that the complete density matrix 
can be determined from the scanned distribution P{x). This is the principle of single scan 
tomography. 

The measurement operator can be derived by projecting the operator n(x) describing 
the position measurement with outcome x into the Hilbert space of the spatial qudit. 

M{x) = ^|si)(si|n(x)|sj)(s,-| = ^M,,(x)|si)(sj|. (2) 

Since x is a continuous variable, all realistic measurements will have a finite resolution. This 
resolution can be described by the measurement operator n(x). In the present experiment, 
the position measurement is realized by a slit of width 6, so that II(x) is given by the integral 

1 /•''/^ 

n(x) = - / \x + x'){x + x'\dx'. (3) 

If (pi{x) = {x\si) is the transverse wave function originating from the i-th slit, Mij{x) can 
be expressed by 

1 /•''/^ 

Mij{x) = - ip*{x + x')(pj{x + x)dx. (4) 

J~b/2 

In the ideal case with infinite resolution, n(x) is the projection operator and the 

idealized measurement operator element is given by 

\im Mij{x) = if*{x)^j{x). (5) 

b—*0 

These idealized patterns have the highest possible visibility. The finite resolution reduces 
the visibility of the measurement. However, quantum state reconstruction only depends on 
the linear independence of the patterns. Therefore the realistic pattern functions of Eq. 
(jl]) result in a more precise reconstruction of the density matrix than the idealized pattern 
functions of Eq. ([5]). 
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FIG. 1: Measurement between the focal and image planes. 



III. MEASUREMENT IN THE INTERMEDIATE PLANE 



The measurement of the spatial qudit is performed by detection in the intermediate plane 
between the focal and image planes of a lens. The lens of focal length / is placed at a distance 
of L from the slit and the photon is detected in the plane at a distance z from the lens, 
where z is between the focal plane and the image plane of the lens. Figure [1] shows the 
schematic setup in the case of qutrits, d = 3, with the slit basis given by {|/), |c), |r)}, each 
element of which corresponds to the photon passing through the left, center, or right slit, 
respectively. 

The wave function ipi{z] x) can be calculated using the Fresnel-Kirchhoff diffraction inte- 
gral for an effective distance of i? = {Lf + zf — Lz) /{z — f) from the multi-slit [l^]. If the 
detection plane is far enough from the image plane, so that the slit width a is sufficiently 
smaller than V RX, where A is the wave length of the photon, the wave function takes the 
form of the conventional slit diffraction pattern given by the sine function sinc(x) = sinx/x, 



ipi{z;x) 



'-Kx 



smc 



K ( x+^-j^r, 



(6) 



Here, K = naf / XR{z — f) scales the sine function and rj is the distance between the optical 
axis and the i-th slit [17|. According to Eq. ([T]), the detection probability is 



P{x) = piiMii{x)+ 2Re [pic] Re [Mi,{x)] + 2Im [p^] Im [Mi,{x)] 

+2Re [pir] Re + 2Im [pi^] Im [M;,(x)] + PccM,,{x) 

+2Re [pcr] Re [M„.(x)] + 2Im [p^r] Im [Mer(x)]+ p„M^r(x). 



(7) 
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In the ideal infinite resolution case for a spatial qutrit, the pattern functions are given by 
Eq. ([5]). In the realistic case of finite resolution, the pattern functions are given by the 
integral in Eq. (jlj) that describes the reduced visibility of the interference patterns. 

To illustrate the effects of finite resolution in the detection of the photon position, Fig. 
[2] shows the pattern functions for infinite and finite resolution of the detector system. The 
dashed lines correspond to ideal measurements with infinite resolution and the solid lines 
correspond to finite resolution measurements. The parameters in the figure were chosen 
to correspond to the experiments described in Section IIVI The detector slit width of 20 
fim that defines the resolution is indicated in the first of the nine graphs. The main effect 
of the finite resolution is to reduce the intensity of the oscillating pattern functions. The 
higher the frequency of oscillation, the more the intensity is reduced. Since the diagonal 
element patterns Ma do not oscillate, their change is negligible. In general, the visibilities 
of interference between slits are weakened more if the slits are farther apart. However, the 
determination of the correct density matrix elements only depends on the use of the correct 
pattern functions, not on their visibilities. 



IV. EXPERIMENTAL SETUP AND VERIFICATION OF ENTANGLEMENT 

In our experiments, we demonstrate the spatial qutrit analysis using entangled photon 
pairs generated by SPDC. A lens images the transverse momenta of the down-converted 
photons onto a pair of triple-slits. After the triple-slits, the photon pairs should ideally be 
in the maximally entangled qutrits states given by 

\^)ab = -^{\l)A\r)B + \c)a\c)b + \r)A\l)B), (8) 

where the suffixes A and B denote the photons in the different paths. 

The experimental setup is shown in Fig. [31 The entangled photon pairs were generated 
by pump beam from a 45 mW cw laser with 405 nm wavelength incident on a 5mm-thick 
/3-barium borate (BBO) crystal. The BBO cristal was set for Type II SPDC in the collinear 
condition. The photons were separated into two different arms by a polarizing beam splitter 
(PBS). The lens before the PBS focuses the momentum eigenstates of the photons on a 
pair of triple-slits, one in each path. The slit width was 45 fim and the distance between 
the slits was 135 fim. To measure the spatial qutrits generated by the triple-slits, a lens of 
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FIG. 2: The nine pattern functions Mij{x) for a spatial qutrit. Dashed hnes correspond to ideal 
measurements with infinite resolution. Solid lines correspond to finite resolution measurements. 
The slit width of the detector system indicated in the upper left hand graph defines the resolution 
of the realistic measurements. 



focal length / = 50 mm was placed at a distance of L = 2/ from each of the triple slits. A 
detector system was placed at a distance of z from each of the lenses. The detector system 
was constructed as shown in inset of Fig. [31 Single slits of width 20 or 40 fim were used 
to select the transverse positions of the photons and an objective lens behind a band pass 
filter of band width 810±5 nm coupled the photons into the multimode fiber connected to 
a photon detector (Perkin-Elmer SPCM-AQR-14) recording the coincidence counts between 
the two arms. Cylindrical lenses were used to reduce the coupling loss between the single 
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slit and the objective lens caused by diffraction at the 20 /im single slits. 

To verify whether our photon source produces the intended entanglement given by Eq. 
([H]), we first confirmed the which-path and interference correlations in the image and focal 
planes of the lenses between the triple-slits and the detector systems. Since the resolution 
was not critical for this purpose, a detector slit width of 40 fim was used in these experiments. 
The which-path correlation was observed in image plane measurements. The detector po- 
sition X in arm A was scanned, while the detector in arm B was fixed at the center of one 
of the slit images. The results are shown in Fig. HI In each scan, the peak appears only at 
the position corresponding to the fixed detector position in arm B, confirming the expected 
correlations between the spatial qutrits. The interference correlation was observed in focal 
plane measurement. In this case, the position in arm B corresponds to a superposition of all 
three slits with a specific phase relation selected by the choice of position. The interference 
patterns obtained are shown in Fig. [51 The phases of the interference patterns clearly de- 
pend on the detector position in arm B, again confirming the expected correlations between 
the spatial qutrits. These two measurements verify that the state generated by our photon 
source in is indeed close to the entangled state given by Eq. (IHl). 



V. DENSITY MATRIX RECONSTRUCTION BY SINGLE SCAN TOMOGRA- 
PHY 



To obtain the complete density matrix of the entangled qutrits, we now apply single 
scan tomography to our photon source. For this purpose, we have to select an intermediate 
plane between the focal and image planes. When the detection plane is near the image 
plane, which-path information corresponding to diagonal elements of the density matrix 
is observed most clearly. Interference information corresponding to off-diagonal elements is 
easier resolved in measurements near the focal plane. In principle, there should be an optimal 
distance z for single scan tomography. Based on our previous work 17|, we started from a 
distance oi z = 1.8/, but found that some improvements in the results of the tomography 
could be achieved by moving the detectors to a distance of z = 1.81/ from the lens. The 
pattern functions shown in Fig. [2] of Sec. IIIII were calculated using this distance between 
lens and detector. 

As the discussion in Sec. Illlt single scan tomography is sensitive to the resolution of 
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FIG. 3: Experimental setup. The transverse momentum is selected by the triple slit in the focal 
plane of the first lens. Photons are separated into two arms by the PBS. The second lens is placed 
at a distance of 2/ from the triple slit in each arm. Coincidence measurements are performed using 
the detector system shown on the right. 
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FIG. 4: Coincidence count rates as a function of scanning detector position for which-path cor- 
relations measured in the image plane. The triangles, dots, and squares show the experimental 
data corresponding to the fixed detector positions at the left, center, and right slit image in arm 
B, respectively. The dashed lines are guides for the eyes. 
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Detector Position ((im) Detector Position (fim) 

FIG. 5: Coincidence count rates as a function of scanning detector position for interference corre- 
lations measured in the focal plane. The dots, squares, and triangles show the experimental data 
and the dashed lines are guides for the eyes, (a) The detector in arm B was set to the positions 
corresponding to phase differences of — vr/S (triangles), (dots), and tt/S (squares), (b) The detec- 
tor in arm B was set to the positions corresponding to phase differences of — vr/G (triangles), vr/G 
(squares), and tt/2 (dots). 



position detection. To achieve as high a resolution as possible, the single slit in the detector 
system was changed to 20 fim and the cylindrical lens was placed behind the slit to com- 
pensate the diffraction effects. However, there remains a non-negligible loss of visibility due 
to the finite slit width, as shown by the pattern functions in Fig. [21 

We have investigated the effects of finite resolution by performing single scan tomography 
of conditional denisity matrices in arm A generated by measurements at fixed detector 
positions in arm B. A comparison of single scan tomography using idealized pattern functions 
neglecting the slit width with tomography using realistic pattern functions including the slit 
width is shown in Fig. O Here, the fixed detector in arm B was near the optical axis at fim 
jisj l. Figure E] (a) shows the coincidence counts obtained by scanning the detector position in 
arm A. We performed two density matrix reconstructions, one using the idealized patterns 
for negligible slit width b^O defined by Eq. ([5]) and the other using the realistic pattern 
functions defined by Eq. (j4]). The dashed line shows the fit with the idealized pattern 
functions and the solid line shows the fit with the realistic pattern functions. Interestingly, 
the fit to the details of the experimental data improves significantly when the finite resolution 
is taken into account. Thus, the application of realistic pattern functions is not only a matter 
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FIG. 6: Comparison of single scan tomography with and without sHt width, (a) Coincidence 
counts obtained in arm A when the fixed detector is near the optical axis in arm B. The dashed 
line shows the fit using the idealized pattern functions for negligible slit width. The solid line shows 
the fit using the realistic pattern functions for the slit width of 20 /im used, (b) Density matrix 
reconstructed with idealized pattern functions, (c) Density matrix reconstructed with realistic 
pattern functions. 



of quantitative errors, but also influences the correct identification of characteristic features 
with elements of the density matrix. This effect is illustrated by the comparison of the 
density matrices reconstructed from idealized pattern functions shown in Fig. Mjo) with 
the more precise density matrix reconstructed from realistic pattern functions shown Fig. 
M^c). The density matrices differ in the values obtained for the off-diagonal elements. In the 
element pir, the change not only affects the size but also the phase of the reconstructed value. 
Thus, the patterns Mir{x) have a rathere subtle effect on the shapes of the reconstructed 
density matrices. 

We performed a sufficient number of conditional density matrix reconstructions for dif- 
ferent setting in arm B to reconstruct the complete 9x9 density matrix of the composite 
system, pab- Figure [7] shows a selection of the conditional scans obtained at different de- 
tector positions in arm B, together with the density matrices reconstructed using realistic 
pattern functions. These scans provide an experimental illustration of the relation between 
the scanned distribution and the density matrix elements. In total, we performed 29 single 
scan tomographies with different detector positions in arm B. We could then determine the 
precise values of experimental parameters such as L by least square fits. The reconstructed 
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density matrix of the complete entangled two qutrits system was found to be 



1 

100 



O.Si 
0.6i 
0.2i 



/ 1.7 -0.1 -O.Si 1.0 + 0.4i 0.5- 

-0.1 + O.Si 4.8 -0.5 + O.Si 1.0- 

1.0 - 0.4i -0.5 - O.Si 28.1 0.2- 

0.5 + O.Si 1.0 + 0.6i 0.2 + 0.2i 4.7 

0.3 + 0.4i 0.7 - l.Oi 29.6 + 5. Si 0.3 - l.Oi 

0.8 + 0.7i 0.1 - O.Si -0.3 - 0.2i -1.4 + 0.7i 

-0.4 - O.Oi -1.0 + O.Oi 24.7 + 4. 7i -0.4 ~ O.Oi 

-1.3 - 0.7i l.l+0.7i -0.5 - O.Si -0.5 + O.Oi 

V 3.2 + 0.2i 0.1+O.li -0.4 + O.li 0.4 + l.Oi 



0.3 - 0.4i 
0.7 + l.Oi 

29.6 - 5. Si 
0.3 + l.Oi 

22.3 
0.6 + 0.4i 

29.7 - 0.7i 
0.7 + 0.3i 
0.2 - 0.4i 



O.S - 0.7i 

0.1 + 0.3 
-0.3 + 0.2: 
-1.4 - 0.7i 

0.6 - 0.4i 
4.1 
-0.3 - 1.5i 
-2.4 + 0.7i 

0.3 - 0.3i 



-0.4 + O.Oi 
-1.0 + O.Oi 

24.7 - 4.7i 
-0.4 + O.Oi 

29.7 + 0.7i 
-0.3 + 1.5i 
27.4 

0.1 - O.Si 

0.6 - 0.9i 



-1.3 + 0.7i 3.2-0.2 
1.1 - 0.7i 0.1 - 0.1 
-0.5 + O.Si -0.4 - 0.1: 
-0.5 - O.Oi 0.4 - 1.0 



0.7 - 0.3i 
-2.4 - 0.7i 
0.1 + O.Si 

5.0 
0.4 - 0.4i 



0.2 - 
0.3 - 
0.6 - 



0.4' 
0.3 
0.9' 



0.4 + 0.4i 
1.9 



Figure [8] shows a comparison of this density matrix with the ideal pure state density matrix 
of the maximally entangled state given by Eq. ([8]). The fidelity of the ideal state \'^)ab was 
= ab{'^\pab\'^)ab = 0.819. Thus the inclusion of position resolution defined by the slit 
width in the pattern functions permits a reliable reconstruction of the quantum state in the 
nine dimensional Hilbert space of the entangled qutrit pair. 



VI. CONCLUSIONS 



Single scan tomography is a method of state reconstruction for an arbitrary dimensional 
spatial qudit that directly connects the scanned distribution to the density matrix elements. 
We can visually understand the elements of the measurement operators as pattern functions 
and a fit of the scanned distribution using these pattern functions provides the reconstructed 
density matrix elements. A reliable reconstruction is possible whenever the pattern functions 
associated with different density matrix elements are linearly independent. It is therefore 
possible to include the effects of finite measurement resolution in the pattern functions 
without affecting the reconstruction of the density matrix. The results obtained with realistic 
pattern functions based on measurement operators describing the finite resolution of the 
experiment are therefore better than the ones obtained with idealized pattern functions 
based on infinitely precise measurement. 

We have demonstrated the effect of including the finite resolution in the pattern function 
on the experimental reconstruction of spatial qutrits. The results show that a reliable 
reconstruction of the density matrix is possible even though the finite spatial resolution of 
20 yum has a non-negligible effect on the pattern functions. In particular, a fidelity of 0.819 
could be obtained for the entangled qutrit pair when the realistic pattern functions were 
used in the reconstruction of the density matrix. 
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FIG. 7: Coincidence counts data. Each measurement was obtained by placing the fixed detector 
in a different position. The sohd hues are fitted to the data points using reahstic pattern functions 
icluding the effects of finite slit width. The corresponding conditional density matrix is shown on 
the right side each interference pattern. 

As the dimensionality of spatial qudit systems increases, it will be increasingly difficult 
to realize high resolution measurements. The method of single scan tomography with finite 
measurement resolution may therefore play a significant role in the realization of higher 
dimensional qudit systems. 
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FIG. 8: Real and imaginary part of (a) the density matrix of the maximally entangled pure state 
and (b) the density matrix reconstructed from single scan tomographies with realistic pattern 
functions. 
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